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Abstract The connection between cosmological observations and neutrino physics is discussed 
in detail. Neutrinos decouple from thermal contact in the early Universe at a temperature of 
order 1 MeV which coincides with the temperature where light element synthesis occurs. Obser- 
vation of light element abundances therefore provides important information on such properties 
as neutrino energy density and chemical potential. Precision observations of the cosmic mi- 
crowave background and large scale structure of galaxies can be used to probe neutrino masses 
with greater precision than current laboratory experiments. In this review I discuss current 
cosmological bounds on neutrino properties, as well as possible bounds from upcoming mea- 
surements. 
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1 Introduction 

In the past few years a new standard model of cosmology has been established 
in which most of the energy density of the Universe is made up of cold dark 
matter and a component with negative pressure, generically referred to as dark 
energy. This model provides an amazingly good fit to all observational data with 
relatively few free parameters and has allowed for stringent constraints on the 
basic cosmological parameters. 

The precision of the data is now at a level where observations of the cosmic 
microwave background (CMB), the large scale structure (LSS) of galaxies, and 
type la supernovae can be used to probe important aspects of particle physics 
such as neutrino properties. Conversely, cosmology is now also at a level where 
unknowns from the particle physics side can significantly bias estimates of cos- 
mological parameters. 

The prime example of this interplay between particle physics and cosmology is 
the use of cosmological data to probe neutrino physics. Particularly the possibility 
to constrain the neutrino mass using cosmological measurements has been studied 
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extensively. 

If only the three known neutrino flavour states, Uf, are considered, they must 
correspond to three mass eigenstates, fj. The two sets are related by 

Uf = Uui, (1) 

where U is the 3x3 neutrino mixing matrix. 

The neutrino oscillation probability in the general 3x3 oscillation case is quite 
complicated, but always depends on Sm^, i.e. on squared mass differences. The 
combination of all currently available data from neutrino oscillation experiments 
suggests two important mass differences in the neutrino mass hierarchy. The 
solar mass difference of Amf2 — 7.5 — 8.7 x 10~^ eV^ {2a) and the atmospheric 
mass difference lAmlgl ~ 1.7 — 2.9 x 10~^ eV^ {2a) (1). In the simplest case 
where neutrino masses are hierarchical these results suggest that mi ~ 0, m2 ~ 
Awgoiar, and nis ~ Amatmospheric- If the hierarchy is inverted one instead finds 
ma ~ 0, m2 ~ Amatmospheric, and mi ~ Amatmospheric- Howcvcr, it is also 
possible that neutrino masses are degenerate, mi ~ m2 ~ ms 3> Amatmospheric- 
Since oscillation probabilities depend only on squared mass differences, Am^, such 
experiments have no sensitivity to the absolute value of neutrino masses, and if 
the masses are degenerate oscillation experiments are not useful for determining 
the absolute mass scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the 
most robust probe of this overall mass scale. Tritium decay measurements by the 
Mainz experiment have been able to put an upper limit on the effective electron 
neutrino mass of rui,^ = {J2i \Uei\'^'mf)^^^ < 2.3 eV (95% conf.) (2) (note that this 
is an incoherent sum so that cancellation due to phases cannot occur). However, 
cosmology at present provides an even better limit, as will be discussed in detail 
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in the next sections. 

Very interestingly there is also a claim of direct detection of neutrinoless double 
beta decay in the Heidelberg-Moscow experiment (3-5). Neutrinoless double 
beta decay is only possible if neutrinos are majorana particles because it requires 
violation of lepton number. The lifetime of the neutrinoless mode of the decay is 
inversely proportional to the square of the effective mass of the neutrino because 
it involves flipping the spin of the internal neutrino. 

The claimed lifetime of ^^Ge measured by the experiment can be translated to 
a neutrino mass. However, there is a large uncertainty in this translation because 
the involved nuclear matrix element is quite difficult to calculate (see for instance 
(6)). Including the estimated matrix element uncertainty the result corresponds 
to an effective neutrino mass in the 0.1 — 0.9 eV range for the parameter niee = 
(note that, contrary to the tritium measurement, this is a coherent 
sum which means that it can be suppressed by phases in the mixing matrix). If 
this result is confirmed then it shows that neutrino masses are almost degenerate. 

Another important question which can be answered by cosmological observa- 
tions is how large the total neutrino energy density is. Apart from the standard 
model prediction of three light neutrinos, such energy density can be either in 
the form of additional, sterile neutrino degrees of freedom, or a non-zero neutrino 
chemical potential. 

In this review I discuss mainly these two questions. Due to the limited space 
available other interesting aspects of neutrino cosmology, such as leptogenesis, 
have been left out. The interested reader is referred to the very thorough review 
(35). 

The paper is divided into sections in the following way: In section 2 I review 
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the thermal evolution of light neutrinos, particularly aspects related to the de- 
coupling of neutrinos around a temperature of 1 MeV. I also discuss Big Bang 
nucleosynthesis and its relation to neutrino physics. 

I section 3 I discuss the current cosmological data available, and section 4 
contains a review of bounds on neutrino properties from this data. Finally, section 
5 contains a discussion. 

2 The thermal history of Hght neutrinos 
2.1 Standard model 

In the standard model neutrinos interact via weak interactions with charged 
leptons, keeping them in equilibrium with the electromagnetic plasma at high 
temperatures. Below T ~ 30 — 40 MeV e+ and e~ are the only relevant particles, 
greatly reducing the number of possible reactions which must be considered. In 
the absence of oscillations neutrino decoupling can be followed via the Boltzmann 
equation for the single particle distribution function (7) 

f-//p| = C„,„ (2) 

where CcoU represents all elastic and inelastic interactions. In the standard model 
all these interactions are 2 <-> 2 interactions in which case the collision integral 
for process i can be written 

1 r (f-p2 d^V^ d?V>4: 



CcolI,i(/l) = J 



2Ei J 2E2{2TTf 2E^{2'Kf 2Ei{2-KY 

X{27rf 6\pi +P2 -P3 +P4)A(/l,/2,/3,/4)-S|M|22_34^., (3) 

where S'|M|^2-^^34,i is the spin-summed and averaged matrix element including the 
symmetry factor S = 1/2 if there are identical particles in initial or final states. 
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The phase-space factor is A(/i, /2, /s, Ja) = fshi^ - /i)(l - /2) - /i/2(l - /3)(1 - 
h). 

The matrix elements for all relevant processes can for instance be found in 
Ref. (8). If Maxwell-Boltzmann statistics is used for all particles, and neutrinos 
are assumed to be in complete scattering equilbrium so that they can be repre- 
sented by a single temperature, then the collision integral can be integrated to 
yield the average annihilation rate for a neutrino 
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r = ^i9l + 9l)T^ (4) 
where 

{sin^ Ow + + sin^ Ow)^ for Ue 
(5) 
sin^ 9w + {-1 + sin^ 6w)'^ for f^,^ 

This rate can then be compared with the Hubble expansion rate 

to find the decoupling temperature from the criterion H = Fl^-^^^^. Prom this 
one finds that Toive) - 2.4 MeV, Td{i'^,t) ^ 3.7 MeV, when = 10.75, as is 
the case in the standard model. 

This means that neutrinos decouple at a temperature which is significantly 
higher than the electron mass. When e+e~ annihilation occurs around T ~ nie/'i, 
the neutrino temperature is unaffected whereas the photon temperature is heated 
by a factor (11/4)V3. The relation T^/T^ = (4/ll)V3 ~ 0.71 holds to a precision 
of roughly one percent. The main correction comes from a slight heating of 
neutrinos by e+e~ annihilation, as well as finite temperature QED effects on the 
photon propagator (8-22). 
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2.2 Big Bang nucleosynthesis and the number of neutrino species 

Shortly after neutrino decoupling the weak interactions which keep neutrons and 
protons in statistical equilibrium freeze out. Again the criterion H = rL_^ 

' ^ — freeze 

can be applied to find that Tfreeze — 0.5^^^ MeV (7). 

Eventually, at a temperature of roughly 0.2 MeV deuterium starts to form, and 
very quickly all free neutrons are processed into ^He. The final helium abundance 
is therefore roughly given by 

'ill III 'In 



Yp 



(7) 

T~0.2 MeV 



1 + Un/rip 

nn/np is determined by its value at freeze out, roughly by the condition that 
nn/'T-p|r=rfreeze ~ e-("^"-"'f)/^fr««-. 

Since the freeze-out temperature is determined by g^, this in turn means that 
5* can be inferred from a measurement of the helium abundance. However, since 
Yp is a function of both fifo/i^ and it is necessary to use other measurements to 
constrain in order to find a bound on 5* . One customary method for doing 
this has been to use measurements of primordial deuterium to infer fi^^^ and from 
that calculate a bound on g^,. Usually such bounds are expressed in terms of the 
equivalent number of neutrino species, = p/pvo, instead of g^. The exact 
value of the bound is quite uncertain because there are different and mutually 
inconsistent measurements of the primordial helium abundance (see for instance 
Ref. (23, 24) for a discussion of this issue). Some of the most recent analyses are 
(23) where a value of 1.7 < < 3.0 (95% C.L.) was found, (25) which found 
-1.14 < AiV^ < 0.73, and (26) which found that = 3.14l°:^5 at 68% C.L. 
The difference in these results can be attributed to different assumptions about 
uncertainties in the primordial helium abundance. It should be noted that all 
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these results are consistent with Ny = 3.04, the standard model result. If a low 
helium abundance is assumed than the BBN bound severely restricts additional 
relativistic energy density beyond the standard model prediction. A full extra 
neutrino degree of freedom is excluded at high significance. However, that is 
not the case when a higher helium abundance is assumed, and the fact that 
the bound is dependent on assumptions should indicate that it is not a level of 
precision where it can be completely trusted. 

Another interesting parameter which can be constrained by the same argument 
is the neutrino chemical potential, = HvjT (27-30). At first sight this looks 
like it is completely equivalent to constraining N^. However, this is not true 
because a chemical potential for electron neutrinos directly influences the n — p 
conversion rate. Furthermore, it is crucial to take neutrino flavour oscillations 
into account when calculating bounds on the neutrino chemical potential. This 
point is discussed in more detail below. 

2.3 Can neutrinos recoup le? 

Standard model neutrinos interact via exchange of W and Z bosons where the 
coupling strength at low energy {E ^ mw,z) is given by Fermi's constant 



The 1/m^ comes from the vector boson propagator. However, at energies beyond 
mw the propagator will have mw replaced by momentum exchange k. Since this 
should be roughly equal to the thermal energy of an average particle in the 
medium, the effective interaction rate should go from 



Gf = 



ira 1 



(8) 



-v/2 sin^ 9wfn^ 



(9) 
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to 



r~ar (r»mvF). 



(10) 



Comparing this with the Hubble parameter which in the radiation dominated 
regime is given by Eq. (6), one finds that 



At high temperatures V/H will be less than 1 and neutrinos will be out of equilib- 
rium. The temperatures where this happens is roug 
which is far beyond the range where the standard model applies. 

However, even if the above example is therefore not physically relevant it shows 
that neutrinos which interact via exchange of massless particles will always have 
an interaction rate (for 2 2 processes) which is 



where is a dimensionless quantity. Therefore such interactions will always drive 
neutrinos toward thermal equilibrium at late times instead of out of equilibrium. 

One example which has received significant attention recently is the possibility 
that neutrinos couple to a massless scalar or pseudoscalar particle. If the inter- 
action is sufficiently strong the neutrinos will come into equilibrium after BBN 
but before the present and consequently they will annihilate and disappear when 
T ~ rui, (31-33). In Fig. 1 we show F/H fov massless standard model neutrinos. 
As can be seen, neutrinos decouple from thermal equilibrium when T ^ 1 MeV. 
However for T S> mz the cross section drops and neutrinos are out of equilibrium 
at very early times. In the same figure we also show what happens if neutrinos 
couple to a light scalar, in which case they come into equilibrium at late times. 




(11) 



T = f3T, 



(12) 
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2.4 The effect of oscillations 

In the previous section the one-particle distribution function, /, was used to 
describe neutrino evolution. However, for neutrinos the mass eigenstates are not 
equivalent to the flavour eigenstates because neutrinos are mixed. Therefore the 
evolution of the neutrino ensemble is not in general described by the three scalar 
functions, /j, but rather by the evolution of the neutrino density matrix, p = ipip^, 
the diagonal elements of which correspond to /j. 

For three-neutrino oscillations the formalism is quite complicated. However, 
the difference in Ami2 and Am23, as well as the fact that sin 26*13 1 means 
that the problem effectively reduces to a 2 x 2 oscillation problem in the standard 
model. The effect of oscillations on the decoupling of neutrinos was discussed in 
(21,22). A detailed account of the physics of neutrino oscillations in the early 
universe is outside the scope of the present review, however an excellent and very 
thorough discussion of this topic can be found in Ref. (35). 

Without oscillations it is possible to compensate a very large chemical potential 
for muon and/or tau neutrinos with a small, negative electron neutrino chemical 
potential (27). However, since neutrinos are almost maximally mixed a chemical 
potential in one flavour can be shared with other flavours, and the end result is 
that during BBN all three flavours have almost equal chemical potential (36-40) . 
This in turn means that the bound on i/g applies to all species. The most recent 
bound on the neutrino asymmetry from BBN comes from (41) 

-0.04 < = M < 0.07 (13) 

for i = e, fi, T. 

The bound assumes complete flavour equilibration during BBN, which with 
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the measured mixing angles and mass differences is a fairly good approximation. 
In models where sterile neutrinos are present even more remarkable oscillation 
phenomena can occur. However, I do not discuss this possibility further, except 
for the possibility of sterile neutrino warm dark matter, and instead refer to the 
review (35). 

2.5 Low reheating temperature and neutrinos 

In most models of inflation the universe enters the normal, radiation dominated 
epoch at a reheating temperature, Trh, which is of order the electroweak scale 
or higher. However, in principle it is possible that this reheating temperature is 
much lower, of order MeV. This possibility has been studied many times in the 
literature, and a very general bound of Trh ^ 1 MeV has been found (42-45) 

This very conservative bound comes from the fact that the light element abun- 
dances produced by big bang nucleosynthesis disagree with observations if the 
universe if matter dominated during BBN. However, a somewhat more stringent 
bound can be obtained by looking at neutrino thermalization during reheating. 
If a scalar particle is responsible for reheating then direct decay to neutrinos is 
suppressed because of the necessary helicity flip. This means that if the reheating 
temperature is too low neutrinos never thermalize. If this is the case then BBN 
predicts the wrong light element abundances. However, even if the heavy parti- 
cle has a significant branching ratio into neutrinos there are problems with BBN. 
The reason is that neutrinos produced in decays are born with energies which 
are much higher than thermal. If the reheating temperature is too low then a 
population of high energy neutrinos will remain and also lead to conflict with 
observed light element abundances. Recent analyses have shown that in general 
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the reheating temperature cannot be below a few MeV (46) (see also (47)). 
3 Cosmological data 

Even though BBN considerations do provide interesting bounds on neutrino 
physics attention has shifted towards observations of the late-time Universe. In 
recent years observations of physics after recombination has provided very strong 
constraint on many cosmological parameters and also on neutrino physics. In 
this section I review the current status of cosmological observations. 

Large Scale Structure (LSS) - The Sloan Digital Sky Survey (SDSS) has 
measured redshifts of close to 1 million galaxies, providing the so far most detailed 
map of the large scale structure of the Universe (48,49). The SDSS data has 
been used extensively for cosmological parameter estimation, including neutrino 
mass bounds. The somewhat smaller 2dFGRS (2 degree Field Galaxy Redshift 
Survey) (50) has also been used for the same purpose. Both surveys provide 
a precise measurement of the galaxy-galaxy power spectrum, Pg{k), down to 
k ~ 0.01/i/Mpc. In order to derive the matter power spectrum, Pm{k), from this 
it is necessary to know the bias parameter, b, defined as 6^(fc) = Pg{k)/Pm{k). In 
general the bias parameter is scale dependent and a function of the complicated 
hydrodynamics of non-linear structure formation. However, if only data on scales 
larger than about k ~ 0.15/i/Mpc is used the bias parameter can be taken to be 
a constant. In principle the bias parameter can be determined from higher-order 
correlations (see for instance (51, 53)), and it does provide a stringent constraint 
on the neutrino mass (52). However the systematic error on the bias parameter 
does rely on the assumption that the non-linear aspects of structure formation 
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are well understood, an assumption which is not necessarily justified at present. 
Therefore the bound on neutrino mass coming from the use of bias constraints 
should probably be regarded as somewhat less robust than bounds relying only 
on the power spectrum shape. 

Cosmic Microwave Background (CMB) — Fluctuations in the CMB were 
first measured by the COBE satellite in 1992 (54). The temperature fiuctuations 
in the CMB can be conveniently decomposed into spherical harmonics 

AT 

— {9,(l>) = J2aimYim{0,(l>)- (14) 

Im 

Prom the a;^ coefficients one can construct the angular power spectrum 

Ct,1 = ihml''), (15) 

where (...) denotes an ensemble average. Since only one realization of the un- 
derlying distribution is available it is in practise replaced by an angular average 
corresponding to averaging over 21 + 1 m-modes for each I. 

The dominant scattering mechanism for photons around the apoch of recombi- 
nation is Thomson scattering on free electrons. Since Thomson scattering polar- 
izes light any anisotropy in the electron distribution leads to a net polarization 
anisotropy in the CMB. Like the temperature anisotropy, polarization can be 
written in terms of angular power spectra. Because there is no circular polar- 
ization component in the CMB there are only two independent components in 
the polarization tensor. The usual decomposition is in terms of curl-free {E) and 
a curl (B) component which yields four independent power spectra: Ct,i, Ce,i, 
Cb,i, and the T-E cross-correlation Cte,i- There is no cross-correlation between 
B and E, T because of different parity. 
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The WMAP experiment has reported data on Ct,i and Cte,i as described in 
Refs. (55-59). Foreground contamination has already been subtracted from their 
pubhshed data. However, it is well known that the WMAP l-year data does 
exhibit some notable anomalies. First, there are several "glitches" in the power 
spectrum, which are presumably due to incomplete foreground removal. Second, 
there is a marked suppression and alignment of the low-w multipoles. There 
has been substantial discussion of the nature of the low-Z anomaly, but with 
no clear conclusion. Very importantly neither of the two anomalies have any 
influence at all on the study of neutrino properties. The reason is that neutrino 
physics affects only scales smaller than roughly the horizon size at recombination, 
and that changes to neutrino physics cannot produce sharp features in angular 
CMB power spectrum. Therefore, any conclusion on neutrinos is likely to hold, 
independent of the possible nature of the anomalies. 

In addition to the WMAP data there are many independent measurements of 
the spectrum on smaller scales by ground based or balloon borne experiments. 
The most important of these at present is the Boomerang experiment (61-63) 
which has measured significantly smaller scales than WMAP. Furthermore, this 
experiment has been the first to measure Ce,i- 

The baryon acoustic peak (BAO) — In the early Universe, prior to recom- 
bination, baryons and photons undergo acoustic oscillations. These oscillations 
are not only detectable in the CMB, but in principle also in the LSS power spec- 
trum. However, since baryons are only a subdominant component early structure 
formation is dominated by cold dark matter, and the oscillations are merely a 
small amplitude modulation of the power spectrum. In terms of the real-space 
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correlation function this modulation corresponds to a well-defined peak, located 
at a scale of roughly lOOh"^ Mpc. This peak was first measured using the bright 
red galaxies of the SDSS (64). The position of the peak provides an accurate 
measurement of the angular diameter distance out to redshift of roughly 0.35. 
This in turn provides a constraint on Qm, ^^DE, and w. Furthermore the con- 
straint shows relatively little sensitivity to other parameters and seems less prone 
to systematics than measurements of the power spectrum amplitude. 

The Lyman-a forest ~ Measurements of the flux power spectrum of the 
Lyman-a forest has been used to measure the matter power spectrum on small 
scales at large redshift. This has been done using high resolution Keck or VLT 
spectra (65, 66). By far the largest sample of spectra comes from the SDSS survey, 
even though the spectral resolution is substantially poorer than in for Keck or 
VLT. In Ref. (67) this data was carefully analyzed and used to constrain the linear 
matter power spectrum. The derived amplitude is A^{k = 0.009 kms~^, z = 3) = 
0.452lQ;Qg and the effective spectral index is n^s = — 2.3211q;q5. Even though 
the higher resolution data extend almost an order of magnitude higher in k the 
smallest scales are likely to be completely dominated by systematics related to 
non-linearities so that relatively little additional information would be gained. In 
order to probe the mass of light neutrinos the SDSS data is therefore optimal. 
On the other hand the small scale data is important for testing models for warm 
dark matter, as will be discussed later. 

The quote SDSS result has been derived using a very elaborate model for 
the local intergalactic medium, including full N-body simulations. It has been 
shown that using the Lyman-a data does strengthen the bound on neutrino 
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mass significantly. However the question remains as to the level of systematic 
uncertainty in the result. Especially the amplitude of the matter power spectrum 
is quite sensitive to model assumptions. Therefore the same caution should be 
applied to bounds using this result as to the bound from LSS data using bias 
constraints, as discussed above. 

Type la supernovae — Measurements of distant type la supernovae provide 
a precise measurement of the luminosity distance di as a function of redshift 
(68,69). This in turn provides an important constraint on H(z) and therefore 
on Cljn, ^DE, and w;de- The most recent and precise data set comes from the 
SuperNova Legacy Survey (SNLS) which has published data for 71 supernovae 
(70) . Another data set which has been extensively used is the Riess et al. " gold" 
data set (71) which consists of 157 supernovae, measured using both HST and 
ground based telescopes. 

Other data — Apart from the above mentioned data the measurement of the 
Hubble constant by the HST Hubble Key Project, i^o = 72 ± 8 km s"^ Mpc~^ 
(72) is important for cosmological parameter analyses. 

4 Neutrino Deirk Matter 

Neutrinos are a source of dark matter in the present day universe simply because 
they contribute to fi^. The present temperature of massless standard model 
neutrinos is T^.^ = 1.95 K = 1.7 x 10~^ eV, and any neutrino with m S> Tj^^o 
behaves like a standard non-relativistic dark matter particle. 

The present contribution to the matter density of Ni, neutrino species with 
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standard weak interactions is given by 



n,h^ = N,^^ (16) 
93.8eV ^ ' 



Just from demanding that Q.,^ <1 one finds the bound (73, 74) 



/46eV 

m^<,—— (17) 



4.1 The Tremaine-Gunn bound 

If neutrinos are the main source of dark matter, then they must also make up 
most of the galactic dark matter. However, neutrinos can only cluster in galaxies 
via energy loss due to gravitational relaxation since they do not suffer inelastic 
collisions. In distribution function language this corresponds to phase mixing of 
the distribution function (75). By using the theorem that the phase-mixed or 
coarse grained distribution function must explicitly take values smaller than the 
maximum of the original distribution function one arrives at the condition 

/CG < /i/,max = -j^ (18) 

Because of this upper bound it is impossible to squeeze neutrino dark matter 
beyond a certain limit (75). For the Milky Way this means that the neutrino 
mass must be larger than roughly 25 eV if neutrinos make up the dark matter. 
For irregular dwarf galaxies this limit increases to 100-300 eV (76, 77), and means 
that standard model neutrinos cannot make up a dominant fraction of the dark 
matter. This bound is generally known as the Tremaine-Gunn bound. 

Note that this phase space argument is a purely classical argument, it is not 
related to the Pauli blocking principle for fermions (although, by using the Pauli 
principle fy < 1 one would arrive at a similar, but slightly weaker limit for 
neutrinos). In fact the Tremaine-Gunn bound works even for bosons if applied 
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in a statistical sense (76), because even though there is no upper bound on the 
fine grained distribution function, only a very small number of particles reside 
at low momenta (unless there is a condensate). Therefore, although the exact 
value of the limit is model dependent, limit applies to any species that was once 
in thermal equilibrium. A notable counterexample is non-thermal axion dark 
matter which is produced directly into a condensate. 

A very interesting direct example of the Tremaine-Gunn bound was studied in 
(78). Here, neutrino clustering in evolving cold dark matter halos was studied 
using the Boltzmann equation. The problem was made tractable by assuming the 
backreaction, i.e. that the CDM halos are not affected by neutrinos. The results 
clearly show that there is a maximum in the coarse grained distribution function 
of 1/2, as expected. The calculation was extended to bosons in (79), where no 
such upper bound was found. In fact bosons can have an average density several 
times higher than fermions with equal mass inside dark matter halos, an effect 
which in principle might be used to distinguish fermionic and bosonic hot dark 
matter. 

4.2 Neutrino hot dark matter 

A much stronger upper bound on the neutrino mass than the one in Eq. (17) can 
be derived by noticing that the thermal history of neutrinos is very different from 
that of a WIMP because the neutrino only becomes non-relativistic very late. 

In an inhomogeneous universe the Boltzmann equation for a collisionless species 
is (80) 

Df _df dx^df dq^df _ 
where r is conformal time, dr = dt/a, and g* = ap^ is comoving momentum. The 
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second term on the right-hand side has to do with the velocity of the distribution 
in a given spatial point and the third term is the cosmological momentum redshift. 

Following Ma and Bertschinger (80) this can be rewritten as an equation for 
^, the perturbed part of / 



f{x\q\T) = fo{q) l + ^{x\q\T 



(20) 



In synchronous gauge that equation is 



Jo OT e amq 



h + 6ri 



-1^ 



Jo 



(21) 



where q^ = qn^ , /x = n^kj, and e = (g^+o^m^)^/^. is the comoving wavevector. 
h and rj are the metric perturbations, defined from the perturbed space-time 
metric in synchronous gauge (80) 



ds^ = a^(r)[-dr^ -|- {Sij + hij)dx'dx^], 



hij = / d-^ke'-" '' ( kikjh{k,T) + {hkj - ■^5ij)Qr]{k,T) 



(22) 
(23) 



Expanding this in Legendre polynomials one arrives at a set of hierarchy equa- 
tions 
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{IFi^i - (/ + l)Fi+i) 



(24) 



For subhorizon scales (/i = r) = 0) this reduces to the form 
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2(7 = 



Fi = 




{lFi_i -{1 + 



(25) 



One should notice the similarity between this set of equations and the evolution 
hierarchy for spherical Bessel functions. Indeed the exact solution to the hierarchy 
is 



This shows that the solution for 6 is an exponentially damped oscillation. On 
small scales, k > r, perturbations are erased. 

This in intuitively understandable in terms of free-streaming. Indeed the Bessel 
function solution comes from the fact that neutrinos are considered massless. In 
the limit of CDM the evolution hierarchy is truncated by the fact that ^ = 0, so 
that the CDM perturbation equation is simply 6 = —h/2. For massless particles 
the free-streaming length is A = cr which is reflected in the solution to the 
Boltzmann hierarchy. Of course the solution only applies when neutrinos are 
strictly massless. Once T ~ m there is a smooth transition to the CDM solution. 
Therefore the final solution can be separated into two parts: 1) A; > t(T = m): 
Neutrino perturbations are exponentially damped 2) A; < t{T = m): Neutrino 
perturbations follow the CDM perturbations. Calculating the free streaming 
wavenumber in a flat CDM cosmology leads to the simple numerical relation 
(applicable only for Teq » m » Tq) (7) 



In Fig. 2 I have plotted transfer functions for various different neutrino masses 
in a fiat ACDM universe (f2^ + Cl,^ + Q,a = 1). The parameters used were 



Fi{kT) ^ ji{kT) 



(26) 




(27) 
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= 0.04, J)cDM = 0.26 - n,,, = 0.7, h = 0.7, and n = 1. 

When measuring fluctuations it is customary to use the power spectrum, P{k, r), 
deflned as 



The power spectrum can be decomposed into a primordial part, Po{k), and a 
transfer function T{k,T), 



The transfer function at a particular time is found by solving the Boltzmann 
equation for (5(r). 

At scales much smaller than the free-streaming scale the present matter power 
spectrum is suppressed roughly by the factor (81) 



as long as Qi, <^ Qm- The numerical factor 8 is derived from a numerical so- 
lution of the Boltzmann equation, but the general structure of the equation is 
simple to understand. At scales smaller than the free-streaming scale the neu- 
trino perturbations are washed out completely, leaving only perturbations in the 
non-relativistic matter (CDM and baryons). Therefore the relative suppression 
of power is proportional to the ratio of neutrino energy density to the overall 
matter density. Clearly the above relation only applies when ili, <^ 0,^, when 
becomes dominant the spectrum suppression becomes exponential as in the 
pure hot dark matter model. This effect is shown for different neutrino masses 
in Fig. 2. 

The effect of massive neutrinos on structure formation only applies to the scales 
below the free-streaming length. For neutrinos with masses of several eV the free- 



P{k,T) = \d\\T). 



(28) 



P{k,T) = Po{k)T{k,T). 



(29) 




(30) 
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streaming scale is smaller than the scales which can be probed using present CMB 
data and therefore the power spectrum suppression can be seen only in large scale 
structure data. On the other hand, neutrinos of sub-eV mass behave almost like 
a relativistic neutrino species for CMB considerations. The main effect of a small 
neutrino mass on the CMB is that it leads to an enhanced early ISW effect. The 
reason is that the ratio of radiation to matter at recombination becomes larger 
because a sub-eV neutrino is still relativistic or semi-relativistic at recombination. 

4.3 Neutrino mass bounds 

4.3.1 Parameter ESTIMATION METHODOLOGY Because massive neutrinos 
affect structure formation it is possible to constrain their mass using a combi- 
nation of cosmological data. The standard approach to cosmological parameter 
estimation is to use Bayesian statistics which provides a very simple method for 
incorporating prior information on parameters from other sources. Using the 
prior probability distribution it is then possible to calculate the posterior dis- 
tribution and from that to derive confidence limits on parameters. There are 
standard packages such as CosmoMC (82), a likelihood calculator based on the 
Markov Chain Monte Carlo method (85,86), available for this purpose. 

It should be noted here that confidence limits derived using Bayesian statistics 
in this way are often much more stringent than implied by frequentist statistics, 
in which one calculates the confidence with which a given hypothesis can be 
excluded. The main reason for this is that observables in our Universe is a single 
realization of an underlying distribution. There is no possibility for producing 
new data to test the hypothesis and therefore the frequentist approach in many 
cases does not fully use the information provided by the given data. Much more 
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detailed discussions of this point can be found in (83,84). In the remainder of 
this section we discuss only results derived using the Bayesian method. 

Independent of the statistical method used results will in general depend on the 
number of parameters used to fit the data. Because of parameter degeneracies 
bounds on a given parameter will in general get weaker if more parameters are 
used in the fit. However, the question remains as to how many parameters should 
plausibly be included. If nothing is known a priori about the cosmological model 
the best-fit for A'^ parameters should in general be substantially lower for A''-!-! 
parameters if the additional parameter is to be included. This point has been 
extensively discussed in (83), and if used at face value indicates that the vanilla 
ACDM model with n = 1 is the preferred model. This model assumes spatial 
flatness and uses only the following free parameters: J7m, the matter density, fi^, 
the baryon density, Hq, the Hubble parameters, r, the optical depth to reion- 
ization, and A, the amplitude of primordial density fluctuations. It furthermore 
assumes that the primordial perturbations are purely adiabatic in nature. 

However, even though the best flt for neutrino mass is for m„ = that 
certainly does not mean the neutrino masses should not be included in parameter 
fltting. The neutrino mass is known from oscillation experiments to be non-zero 
and therefore must be included in parameter estimation. 

Furthermore almost all studies of constraints on non-standard cosmological 
parameters assume the minimal ACDM model as the benchmark and then in- 
clude the non-standard parameters specific to the given study. This approach is 
dangerous because there could easily be several non-standard parameters which 
produce a much better fit when combined. It is important to study more gen- 
eral models, especially since numerical techniques have now made it feasible to 
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use many more free parameters. Below we will discuss one case where including 
two non-standard parameters can lead to completely changed conclusions about 
neutrino mass bounds. 
4.3.2 Current bounds ^ 

With the WMAP data alone sensitivity to the neutrino mass is relatively lim- 
ited. In a standard likelihood analysis for the 6 parameters of the minimal ACDM 
model plus neutrino mass we find an upper bound of J2 i^u < 2-1 eV (95% C.L.). 
This value is consistent with that derived in other recent analyses (87, 88) . 

Including LSS data the bound can be improved. However, if there is no in- 
formation on the normalization of the large scale structure power spectrum (the 
bias parameter), the improvement is relatively modest. For the same model as 
above the bound for WMAP+SDSS is 1.8 eV (49). 

Once other data, such as type la supernova measurements, are included the 
bound can be improved significantly. A very stringent constraint comes from 
including information on the normalization of the matter power spectrum from 
measurements of the Lyman-a forest. In (52) a bound of X^mj, < 0.42 eV (95% 
C.L.) was derived using WMAP, SDSS, SNI-a, and Lyman-a data for the stan- 
dard ACDM model with neutrino mass ((89) found an almost identical result). 
However, the uncertainty related to the modelling of the intergalactic medium at 
high redshift does make the bound model dependent. 

Another very interesting bound can be derived when the measurement of the 

baryon acoustic oscillation (BAG) peak (64) is included (91). Even for a very 

general model with 11 free parameters (adding the number of neutrino species 

N^, a possible running of the spectral index, a^, and the dark energy equation 
Another recent review discussing neutrino mass bounds is (90) 
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of state, w) the upper bound was found to be 0.48 eV (95% C.L.) (91). This is 
quite interesting because without the BAO data it was shown in (92) that there 
is a severe degeneracy between X^mjy and w. For CMB, LSS and SNI-a data the 
bound on X^m,/ can be relaxed by as much as a factor of 3 when w is allowed 
to vary. This degeneracy is broken by adding the BAO data because it provides 
a tight relation between and w. The result is also encouraging because the 
systematic uncertainty in the BAO measurement is likely to be much smaller than 
in the Lyman-a forest data. Furthermore the derived best fit model is almost 
exactly identical to the best fit model found in (52) , indicating that the Lyman-a 
and BAO data are completely consistent. 

Finally, if both BAO and Lyman-a data is combined the bound strengthens to 
J2 < 0.30 eV (95% C.L.) (91), but with the same caveat with regards to the 
Lyman-a data as before. 

In conclusion the current bound on the sum of neutrino masses can be in the 
range between 0.3 and 2 eV, depending on the data and parameters used. In 
general: 

a) The bound on for the minimal ACDM plus neutrino mass (8 parameters 
in total) is in the 0.5-0.6 eV range if CMB, LSS, and SNI-a data is used. 

b) This bound can be relaxed by a large factor when more parameters, such as 
w and Ny are included, showing that the bound is not robust. 

c) Even including several new parameters the bound can be pushed below 0.5 eV 
when additional data is added to the basic CMB, LSS, and SNI-a data. Either 
baryon acoustic oscillation or Lyman-a forest data can be used, with almost 
identical result. Especially the BAO data seem to be almost free of systematics. 
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and consequently the bound seems to be robustly in the 0.5 eV range. 

For those interested an (incomplete) list of other studies of neutrino mass bounds 
from cosmological data can be found in (93-98). 

In Table 1 various upper bounds on neutrino mass are summarized. The cases 
1-4 are taken from (91) and show quite clearly that results are strongly dependent 
on data and parameters used. Also shown in this table is a selection of other 
recent results, all based on the standard ACDM plus neutrino parameter set. 

In Fig. 3 we show the upper bound on the sum of neutrino masses for various 
cases which are presented in Table 1. Also shown in this figure is the sum of 
neutrino masses as a function of the mass of the lightest mass eigenstate. The 
red band is for the normal hierarchy and the blue is for the inverted hierarchy. 
All bounds were calculated assuming that the neutrino mass is shared equally 
between all species. As can be seen from the figure that assumption is clearly 
justified at present. In the same figure we also show the best fit region for the 
claimed detection of neutrinoless double beta decay by the Heidelberg-Moscow 
experiment (3-5). There does seem to be some tension between the derived upper 
bound from cosmology and this result. However, at present it seems premature to 
exclude it based on cosmological arguments. This is particularly true because the 
exact value of the mass inferred from the experiment is highly uncertain because 
of the uncertainty in the nuclear matrix element calculation. 

In this section we have discussed only active neutrino species. However, the 
framework applies equally well to sterile neutrinos and other light species. For 
instance (99) have derived stringent limits on light sterile neutrinos (see also 
(100)). 
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4.3.3 Neutrino interactions beyond the standard model The mass 
bounds derived from cosmology depend on the assumption that neutrinos have 
no interactions beyond the standard model weak interactions. However, exper- 
imental constraints on non-standard neutrino interactions are in general quite 
weak and allows for the possibility of new, cosmologically important interactions. 
One example is neutrinos coupling to a light scalar or pseudoscalar, as was dis- 
cussed in Section 2.3. In this case neutrinos come into thermal equilibrium with 
the scalar late in the evolution of the Universe. This was studied in (31) as a 
possible way of circumventing the cosmological bound on neutrino masses. The 
idea is that neutrinos coupling to a new massless degree of freedom will anni- 
hilate and disappear once T ~ rrii,, leaving only little imprint on large scale 
structure. However, neutrinos which are strongly coupled will undergo acous- 
tic oscillations prior to CMB formation instead of damping via free streaming 
(anisotropic stress). Since neutrino perturbations act as a source term for the 
baryon-photon oscillations this will increase the amplitude of CMB anisotropies 
on scales small than the horizon scale at recombination. This effect was studied 
numerically in (32, 101, 102). Both studies find that this effect disfavours strongly 
coupled neutrinos, although the exact level of significance differs. The same ef- 
fect was also studied analytically in (103-105). It should be noted that (32, 102) 
studied the effect only in the strong coupling limit. In the intermediate regime 
the situation is more comphcated (see (33)). 

The fact that strongly interacting neutrinos are disfavoured by CMB and LSS 
data means that any model producing neutrino interactions strong enough to 
fast momentum exchange around the epoch of recombination (34). In the case 
of interactions with a light pseudoscalar like the majoron this leads to a very 
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rough bound on the diagonal couphng ga of 10~^ from consideration of pair 
creation and annihilation processes. For the off-diagonal coupling it is roughly 
gij < 10~^^(0.05eV/m)^ which is in fact the strongest known bound on these 
parameters. An exact calculation of the bound would require a solution of the 
Boltzmann equation for neutrinos in the intermediate coupling regime around 
recombination. 

Another possibility which has received quite a lot of attention recently is the 
possibility of neutrinos coupling to a very light scalar field (106) (see also (107, 
108)). In this case the coupling can produce an effective neutrino mass which 
is density dependent (mass varying neutrinos), and by tuning the coupling it is 
possible to have the combined neutrino-scalar fluid have an effective equation of 
state close to w = —1, i.e. dark energy. It was subsequently argued that the 
particular form of the potential in (106) leads to a small scale instability which 
drives w to (109) (see also (110)). However, this instability apparently requires 
two conditions to be fulfilled: a) The effective mass of the light scalar should 
have m » H and b) The equation of state should obey u; > — 1. Both these 
conditions can be violated by other potentials (see for instance (111, 112)). Mass 
varying neutrinos have been discussed in other contexts, for instance in (113-116). 
4.3.4 The number of neutrino species - joint CMB and BBN analysis 
The BBN bound on the number of neutrino species presented in the previous 
section can be complemented by a similar bound from observations of the CMB 
and large scale structure. The CMB depends on Ni, mainly because of the early 
Integrated Sachs Wolfe effect which increases fluctuation power at scales slightly 
larger than the first acoustic peak. The large scale structure spectrum depends 
on Ni, because the scale of matter-radiation equality is changed by varying Ni,. 
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Many recent papers have used CMB, LSS, and SNI-a data to calculate bounds 
bounds on N^^ (23, 117-120), and some of the bounds are listed in Table 4.3.4. 
Recent analyses combining BBN, CMB, and large scale structure data can be 
found in (23, 118), and these results are also listed in Table 4.3.4. 

Common for all the bounds is that = Q is ruled out by both BBN and 
CMB/LSS. This has the important consequence that the cosmological neutrino 
background has been positively detected, not only during the BBN epoch, but 
also much later, during structure formation. 

The most recent bound which uses the SNI-a "gold" data set, as well as the 
new Boomerang CMB data finds a limit of = 4.21^;^ at 95% C.L. (121). The 
bound from late-time observations is now as good as that from BBN, and the 
two derived value are mutually consistent given the systematic uncertainties in 
the primordial helium abundance. 

In Fig. 4 we show the currently allowed region for flj^h? and Ni, (taken from 
(121)), showing the overlap between BBN constraints and the CMB+LSS+SNI-a 
constraint. 

In principle, even if the neutrino mass is very small, it could be possible to 
detect the difference between additional relativistic energy density, parameterized 
by A'^i^, and a non-thermal distribution of the active neutrinos (122). However, 
it is unlikely that the precision of observational data from planned experiments 
will be sufficient to measure the difference. 

4.3.5 Future neutrino mass measurements The present bound on the 
sum of neutrino masses is still much larger than the mass difference, |Am23| ~ 
0.05 eV (123, 124), measured by atmospheric neutrino observatories and K2K . 
This means that if the sum of neutrino masses is anywhere close to saturating 
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the bound then neutrino masses must the almost degenerate. The question is 
whether in the future it will be possible to measure masses which are of the order 
Am23, i-e. whether it can determined if neutrino masses are hierarchical. 

By combining future CMB data from the Planck satellite with a galaxy survey 
like the SDSS it has been estimated that neutrino masses as low as about 0.1 eV 
can be detected (125, ?) when analyzed within the minimal 8 parameter ACDM 
model. However, the degeneracy between J^i^v found in present data 

is likely to persist even with much more accurate future data (92), limiting the 
precision with which the neutrino mass can be measured. 

One very powerful probe which will become available in the future is weak 
gravitational lensing. With future CMB experiments like Planck (130) and the 
proposed Inflation Probe (131) it will be possible to measure lensing effects on 
the CMB itself. In this case it seems likely that a sensitivity below 0.1 eV can 
also be reached with CMB alone (127-129). 

Another possibility is to probe weak lensing at much lower redshifts. This can 
be done by measuring lensing of background galaxies by large scale structure. 
This has the major advantage that it is possible to perform lensing tomography 
by binning data in different redshift bins. In (132) the projected sensitivity 
from combining future CMB data with a large scale weak lensing survey was 
studied. It was estimated that at la the sensitivity could be pushed to 0.03 eV. 
Such lensing surveys are planned with experiments such as the Large Synaptic 
Survey Telescope (LSST) which is planned to start in 2012 (135). At present the 
systematic uncertainties related to such surveys is to some extent unknown. One 
source of error is related to the uncertainty in measurements of galaxy shapes 
and another is related to the fact the redshift of galaxies must be estimated 
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photometrically since spectroscopy is not feasible for the large number of galaxies 
needed. Estimates in general show that it will be possible to control the errors 
at the required level (133, 134). 

Another promising tool which could further increase sensitivity is to use future 
cluster surveys (136). Such a survey could also bring the sensitivity to 0.03 eV. 

As noted in Ref. (?) the exact value of the sensitivity at this level depends 
both on whether the hierarchy is normal or inverted, and the exact value of the 
mass splittings. At this level, neutrino masses cannot be regarded as degenerate 
and the normal and inverted hierarchy models must be tested separately. 

4.4 Neutrino warm dark matter 

While CDM is defined as consisting of non-interacting particles which have essen- 
tially no free-streaming on any astronomically relevant scale, and HDM is defined 
by consisting of particles which become non-relativistic around matter radiation 
equality or later, warm dark matter is an intermediate. One of the simplest pro- 
duction mechanisms for warm dark matter is active-sterile neutrino oscillations 
in the early universe (138-142). 

One possible benefit of warm dark matter is that it does have some free- 
streaming so that structure formation is suppressed on very small scales. This 
has been proposed as an explanation for the apparent discrepancy between ob- 
servations of galaxies and numerical CDM structure formation simulations. In 
general simulations produce galaxy halos which have very steep inner density 
profiles p oc r", where a ~ 1 — 1.5, and numerous subhalos (143, 144). Neither 
of these characteristics are seen in observations and the explanation for this dis- 
crepancy remains an open question. If dark matter is warm instead of cold, with 
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a free-streaming scale comparable to the size of a typical galaxy subhalo then 
the amount of substructure is suppressed, and possibly the central density profile 
is also flattened (145-151) . In both cases the mass of the dark matter particle 
should be around 0.5 - 1 keV (153-155), assuming that it is thermally produced 
in the early universe. 

On the other hand, from measurements of the Lyman-a forest flux power spec- 
trum it has been possible to reconstruct the matter power spectrum on relatively 
small scales at high redshift. This spectrum does not show any evidence for sup- 
pression at sub-galaxy scales and has been used to put a lower bound on the 
mass of warm dark matter particles of roughly 500 eV (156). An even more se- 
vere problem lies in the fact that star formation occurs relatively late in warm 
dark matter models because small scale structure is suppressed. This may be in 
conflict with the low-/ CMB temperature-polarization cross correlation measure- 
ment by WMAP which indicates very early reionization and therefore also early 
star formation. One recent investigation of this found warm dark matter to be 
inconsistent with WMAP for masses as high as 10 keV (145). 

Very interestingly a keV sterile neutrino will decay radiatively and produce x- 
ray photons which can reionize the Universe at high redshift. The decay lifetime 
is given purely in terms of the mass and mixing angle with active species 



Since the total contribution to the energy density is also given purely in terms of 
these parameters (assuming a zero lepton assymetry) it is possible to calculate 
the total x-ray photon intensity as a function of redshift. 

Such decays have been proposed as a possible explanation for the very high 
optical depth indicated by the WMAP polarization measurement which seems 




(31) 
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to require partial reionization already at 2; ~ 15 — 20 (157-160). However, in 
this case the mass must be several keV in order to have a short enough lifetime 
(since r oc m~^). Therefore there does not seem to be a common mass which can 
address both the structure formation problems and the reionization problems. 

The case for warm dark matter seems quite marginal, although at present it is 
not definitively ruled out by any observations. 

5 Discussion 

In the present paper I have discussed how cosmological observations can be used 
for probing fundamental properties of neutrinos which are not easily accessible in 
lab experiments. Particularly the measurement of absolute neutrino masses from 
CMB and large scale structure data has received significant attention over the 
past few years. Prom cosmological observations it has been possible to derive an 
upper bound on the sum of neutrino masses which seems to be robustly in the 
0.5 eV range. This is substantially better than the present bound from tritium 
decay measurements which is m,^^ < 2.3 eV (95% C.L.), leading to J2 ^ 7 eV. 

In the future this type of measurement will be improved by more than an 
order of magnitude by the KATRIN experiment (161, 162) which has a projected 
sensitivity of 0.2 eV on the effective electron neutrino mass. 

Another cornerstone of neutrino cosmology is the measurement of the total en- 
ergy density in non-electromagnetically interacting particles. Por many years Big 
Bang nucleosynthesis was the only probe of relativistic energy density, but with 
the advent of precision CMB and LSS data it has been possible to complement 
the BBN measurement. At present the cosmic neutrino background is seen in 
both BBN, CMB and LSS data at high significance, with N,y = being excluded 
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at more than 5a. 

Finally, cosmology can also be used to probe the possibility of neutrino warm 
dark matter, which could be produced by active-sterile neutrino oscillations. 

In the coming years the steady stream of new observational data will continue, 
and the cosmological bounds on neutrino will improve accordingly. For instance, 
it has been estimated that with data from the upcoming Planck satellite it could 
be possible to measure neutrino masses as low as 0.05 eV which would allow for 
a determination of the neutrino mass even in the normal hierarchy. 
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T (eV) 



Figure 1: The interaction rate for standard a standard model neutrino is shown 
in red. The blue curve is for a 0.003 eV neutrino coupled to a massless scalar 
with dimensionless coupling 5 = 3 x 10^^. 
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Figure 2: The transfer function T(A;, t = to) for various different neutrino masses. 
The sohd (black) hne is for ruu = 0, the long-dashed for rrii, = 0.3 eV, and the 
dashed for = 1 eV. 
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Figure 3: The 95% upper bounds on J2 "^^u for the different cases presented in 
Table 1. The horizontal axis is the mass of the lightest neutrino mass eigenstate. 
The set of full red lines correspond to the normal neutrino mass hierarchy and the 
blue to the inverted mass hierarchy. The horizontal yellow band is the claimed 
result from the Heidelberg-Moscow experiment. Note that the exact allowed 
region from this experiment is somewhat uncertain due to uncertainties in the 
nuclear matrix element. Figure taken from (91). 
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Case m„ (95% C.L.) 

1: 11 parameters, CMB, LSS, SNI-a 2.3 eV 

2: 11 parameters, CMB, LSS, SNI-a, BAG 0.48 eV 

3: 8 parameters, CMB, LSS, SNLa, BAG 0.44 eV 

4: 8 parameters, CMB, LSS, SNLa, BAG, Ly-a 0.30 eV 
Other recent bounds 

8 parameters, WMAP only (87) 2 eV 

8 parameters, WMAP, SDSS (49) 1.8 eV 

8 parameters, WMAP, SDSS, SNLa, Ly-a (52) 0.42 eV 



Table 1: Mass bounds using different sets of data and parameters (taken from 
(91). 



Table 2: Various recent limits on the effective number of neutrino species, as well 
as the data used. 



Ref. 


Bound on (95% C.L.) 


Data used 


Crotty et al. (119) 


1.4 < 


< 6.8 


CMB, LSS 


Hannestad (118) 


0.9 < 


< 7.0 


CMB, LSS 


Pierpaoh (120) 


1.9 < N„ 


< 6.62 


CMB, LSS 


Barger et al. (23) 


0.9 < A^^ 


< 8.3 


CMB 


Hannestad (121) 


3.0 < 


< 5.9 


CMB, LSS, SNLa 
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Figure 4: The 68% (dark) and 95% (light) hkelihood contours for ^hh^ and N,y 
for all available data. The other contours are 68% and 95% regions for BBN, 
assuming the ^He and D values given in (26). Figure taken from (121). 



